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Abstract 

Collisionless energy absorption in ID nanoplasma layer is considered. 
Straightforward classical calculation of the absorption rate in action-angle 
variables is presented. In regular regime the result obtained is the same 
as in but deeper insight is possible now due to the technique used. 
Chirikov criterion of the chaotic absorption regime is written out. Col- 
lisionless energy absorption rate in nanoplasma layer is calculated in 
stochastic regime. 
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1 Collisionless absorption in regular regime. 

One of the novel problems of laser-matter interaction is the problem of energy 
absorption in nanometer targets subjected to an ultrashort (up to few picosec- 
ond) intence (10^^ — lO^^W/cni ) laser fields. During such interaction hot (up to 
several keV energy) classical plasma bounded in nanoscale volume is produced, 
which has a life time of about hundreds of femtoseconds. This is dense plasma 
with the electron density of 10^'^cm^^ and more. Such systems are used to be 
called nanoplasma since first experiments of intense short laser interaction with 
three-dimensional nanobodies (atomic Van-der-Vaals clusters) were hold in 1996 

ID 

Nanobodies are known to absorb much more compared to traditional tar- 
gets like gas or even bulk. The great amount of energy contained in tiny volume 
results in breakdown, birth of energetic particles and high harmonics genera- 
tion [21 El ■ Different mechanisms of absorption were suggested to explain such 
phenomena. They are inner ionization, inverse bremsstrahlung effect, vacuum 
heating, collisionless heating and some others a bit more sophisticated. 

As far as nanoplasma is a strongly bounded system with the width much less 
than laser wave length, the most interesting mechanism of energy absorption in 
it is collisionless heating in self-consistent potential. It was considered recently 
in one-dimensional systems corresponded to irradiated films and more deeply in 
three-dimesional systems which correspond to nanoclusters; the important role 
of it in the absorption process was evidently shown (for ID situation see 0]). 
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The problem of coUisionless energy absorption in thin films irradiated by 
intence short laser pulse was considered in ^ in the frames of the following 
model. First, the incompressible liquid approximation for the electronic cloud 
was used^ and both the self-consistent potential and distribution function were 
taken as if they are known function slowly changing on times of the laser pulse 
duration. This means that the self-consisted system of Boltsman and Laplace 
equations was supposed to be solved elsewhere. Then, in |^ dipole aproxima- 
tion^ and the perturbation theory on the small parameter of the dimensionless 
oscillation amplitude of the electron cloud^ was used. Here we use the same 
model. The calculation of the rate of coUisionless absorption presented in 0] 
was based on quantum-mechanical approach in quasiclassical limit. Note, that 
the system considered is classical, and the final result in ^ does not contain 
Plank constant. Althogh this method is non-contradictory, it hides some clas- 
sical features of the system. The present paper has the aim to make the same 
calculation for this system in the frames of classical mechanics and to learn more 
about it possible behaviour. 

Let the particle is bounded in self-consistent potential U{z), with the energy 
distribution function F(e). Mean absorbed energy is defined as: 

d6F(e)g(e), (1) 

where q{e) - is the work of the field over the particle with energy e in the time 
unit, averaged with the initial condition for this particle. Due to the shielding 
effect external 

E(i) = Eocos(wi). (2) 

and internal £ laser fields differs, and in the model of incompressible fluid used 
in |1| there is a relation 

= , (3) 

J(l-w2/^2)2+4r2/^^2 

where ujp is plasma frequency, F is the damphing constant, and we direct electric 
field along z-axis. 

One-particle Hamilton function of the system considered has the form: 

2 

H{p, z, t) = + U{z) - e£^zcos{Lut + a) = Ho + Mnt. (4) 

It is convenient to come to the action- angle variables {1,0), defined for the 
non-perturbed system in a standart way 0: 

dSg{z,I) 



I=^f v/2me(6 - L/(z))dz, 0= ' , (5) 



^For more details about the applicability of this assumption see (JJ. 

^System has linear dimensions much less than the laser wave length, A = SOOnm for typical 
Ti;Sa laser. 

''For the reasonable parameters of the system, such as electronic density lO^^cm""', laser 
intensity 5 ■ lO^^W/sm^, ocsillation amplitude has the order of Inm, and the typical width of 
the film is lOOnm. See W for details. 
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where 

5g(z,Z) = / V2™c(e(I) - U{z))dz (6) 



is a generating function. Using the expansion of the coordinate z{t) in Fourier 
series 

z{I, 0) — ^ Zn{I) cos 710, 
n 

SO that z{I, 8) coswt — ^ z„(Z)[cos(n6 + wt) + cos(n8 — ut)] (7) 

n 

we omit highly oscillating terms. Then Hamilton function looks like 

=Ho(I) - ^^^.(I)cos((2s + l)0-wt), zn = z,. (8) 

s 

Slow change of the one-particle distribution function on large time scale is 
a diffusion in action (energy) space (see, for example, 0), and may be defined 
by the Fokker-Plank-Kolmogorov equation 

gF(6,t) _ld , 9F(6,t) 
dt ~2de^'' de ■ 
Here we first introduce the diffusion coefficient 

D = ^^^, (10) 

which is defined on large (ideally infinite) observation time T, and do not de- 
pend on it. 5T is the addition to unperturbed action under the influence of 
the perturbation in Q. To define it we should write down and solve motion 
equations: 

i = ^{'^s + l)zs{T) sintZ/ 

S 

= n{i) -%^^n{i)z'^[i)cos^, I' ^ {2s + 1)0 - ut. (11) 



Then 



5Ii = ^-^Y.^2s + l)zM-^^^^^^, Ps = {2s + l)n{e,)-u, (12) 

and finally for the diffusion coefficient H10|l we get 

e2^2 ,^, ((1-cos/3,T)(1-cos/3.T))t 
D = -^^^{2s + l){2r + l)z,{l)zr{T) —q ■ 

s r HSHT 

(13) 

To get the energy gain in regular regime without stochasticisy in such description 
we should take the formal limit T ^ oo, because nothing can change the particle 
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trajectory in coUisionless system. During the averaging procedure delta-function 
and delta-symbol Srs appears, and finaly diffusion coefficient in regular regime 
reads as 

D(J) ^ ^ ^(2. + irzUl)nS{Ps). (14) 

s 

To obtain energy gain we integrate FPK-equation @ by parts and come from 
action to energy we get: 



d{e) 1 /^om ^ nm^I^^T W^rn 



This result repeats the result from ^ , where it was obtained on the quasiclassical 
language. The sum over s here means that only the particle at resonant levels 

w = (2s + l)f7(e) (16) 

can absorb energy. 



2 Chirikov criterion of the stochasticity. 

The technic presented in the previous section allows to describe classical motion 
of the particles in nanoplasma in the region of parameters, where quasiclassical 
description used in p| fails. Supposing that the particle has the initial energy 
close to resonance energy defined by the condition (|16ll , averaging over time I0J 
we obtain so-called resonant Hamilton function 

■H{I,e)^n^{I)-^Zs{l)cos{{2s + l)0-ujt), z^ = zs. (17) 

According to let us carry out one more canonical transformation with the 
help of generating function 

g(X,.Z',<) = -(Z-J.)|^, (18) 

with Ts - resonant action, corresponded to the resonance energy with the num- 
ber s from 116|l . For new (P, if') and old (I, 0) variables the following relations 
are fulfilled 

p = -^^^ e = -^ = ^^ (19) 

dW 2s + l' dl 2s + 1 ^ ' 

In new variables Hamiltonian has the form 

n{P, f ) = i(2s + \f^l[{es)n{es)P^ - ^zs{es) cos^, (20) 
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where the decomposition on the small deviation from the resonance action is 
presented. Hamiltonian function H20I) describes the nonlinear mathematical pen- 
dulum. Chaotic motion begins when the amplitude in action of such a system 
is greater than the distance between neighbour resonances^ 0: 




Figure 1 : The resonances overlapping. 



\n'^{es)\n{es) 



>|I(e.)-J(e,+i)|. 



(21) 



This condition means that the action oscillation near one resonance come to the 
region, occupied by the nearest another one (see Fig^. The distance between 
resonances in terms of frequences, according to condition is 



2s + 1 2s + 3 



to 

2?' 



(22) 



Then, with H21() . the criterion of the stochastity appearance, expressed through 
the field strength inside the system is: 

(23) 
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This is overestimation: when the inequality is fulfilled, chaos has to take place. 
In reality threshold is lower [7]. For model potentials from (|23|l one can find 
that in rectangular well of the 100 nm width chaotic regime begins if the inner 
field © is about lOi^W/cm^ 

■^The Chirikov criterion of the resonances overlapping. 
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3 Collisionless absorption in stochastic regime. 



In the situation when chaotic behaviour is developed, to calculate the heating 
rate in H13() it is nesessary to use for averaging not infinite time, but such a 
time Tc which define the dynamical memory of the particle about its previous 
history. The standart procedure to find it is to define mapping of angle variable 
which in our system is 



Qn+i = On + ^^^^ +icyxs sin((2s + l)Qn + as), 



(24) 



where 



(2s + l)7rO 
as — arccos 



Bs 



As 



dzs{I) z^(2s + l) do. 
uj-{2s + 1)VL VV dl ^ (2s + l)f7aX 
Zs{2s + l)7rcos((2s + l)7rf]/w) dVl 



sin((2s + l)7rf7/w) 



Lo dX ^ 

37r(2s + l)zs sin((2s + l)7rf^/tj) dn{I) 



LU - (2s + 1)0 



dJ 



(25) 



For the step time of mapping (|24|l the period of external laser field was taken. 
For such mapping the decorellation time Tc can be estimated as |S] 



47r 



As djlnAT 



(26) 



where As is the number of essential items Xs in sum (|24|l . In our situation it 
has the order of 2s + 1. Finally, the diffusion coefficient can be obtained from 
H13|l with substitution T Tc- It reads 



D(J) = ^Y. II(2s + l)(2r + 1)z,(J)z,(J)A,,(t,), 



(27) 



where 



Arsir) = 



(1 — cos (3st){1 — cos PrT) 

TPsf3r ■ 



(28) 



Energy gain can be obtained from the FPK-equation in the same way as we did 
it earlier for regular regime H15|) : 

^ = D(J) = ^ / d6f7(/(e))^^^^5](2s+l)(2r+l)z,(X(6))z.(J(e))A.,(e,r,). 



dt 



(29) 
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In this expression all resonance levels take part in the absorption process simul- 
taneously. Moreover, in such a situation particle with arbitrary energy should 
gain energy from the external field. Formula H29|l is the main result of the 
present work. It describes the collisionless heating in ID classical nanoplasma 
layer when the field strength is enough for chaotic regime to take place, accord- 
ing to 1(23) • 

Author wishes to thank S.V. Popruzhenko, D.F. Zaretsky, I.Yu. Kostyukov 
for fruitful discussions. The work was done with the financial support of RFBR. 
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